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VISCOELASTIC  ANALYSIS  OF  STRESSES  AND  STRAINS  IN  SOLID 
PROPELLANT  GRAINS  UNDER  PRESSURE  LOADINGS 


Wang  Yuanyou 
Abstract 


This  paper  is  a  viscoelastic  analysis  of  the  stresses  and 
strains  caused  by  internal  pressure  loadings.  Some  engineer¬ 
ing  analysis  methods  already  exist  for  solving  this  classical 
problem.  For  example,  formulas  were  derived  on  the  assumption 
that  propellants  are  incompressible  materials  or  overlooking 
the  casing  deformation.  Although  other  calculation  formulas 
do  not  make  these  assumptions  yet  Schapery's  approximate 
inversions  were  used  to  obtain  viscoelastic  solutions.  This 
paper  proposes  a  new  calculation  method  which  does  not  use  the 
above  assumptions  or  the  approximate  inversions.  Therefore, 
the  solutions  of  the  stresses  and  strains  are  more  accurate 
than  the  above  mentioned  methods.  This  paper  introduces  two 
planar  problems  and  their  results  under  suddenly-applied 
pressure  and  gradually-applied  pressure. 


Symbols 


a 

b 

c 

E 

k 

k. 

1 

P. 


inside  radius  of  grain; 

outside  radius  of  grain; 

outside  radius  of  casing ; 

elastic  constant; 

stress  concentration  coefficient; 

stress  concentration  coefficient  of  inner  hold 

internal  pressure; 


E 

g 

E 

c 

E(t) 
F  ( t) 
F 

H  ( t) 


propellant's 

propellant's 

propellant' s 
propellant 's 
propellant' s 
propellant's 

(t)  ,  t  7/  0 

,  t  <  0, 


l 


glass  modulus; 

rubber  modulus; 

relaxation  modulus; 
creep  compliance; 
rubber  compliance; 
glass  compliance; 

unit  step  function; 


h  :  thickness  of  casing  wall; 

P'  :  contact  pressure; 


r  :  radius  of  grain; 

u  :  radial  displacement 

&  :  (  /l  ,v)  correlation  parameter; 

s  f fi  K  /l  #v,vc) correlation  parameter; 

v  :  propellant's  Poisson  ratio; 
vc  :  Poisson  ratio  of  casing  material. 

7P  :  b/ as  . 

I.  Introduction 

During  the  ignition  of  a  solid  rocket  engine,  internal 
pressure  can  cause  the  inner  surface  of  propellant  grain  to 
have  cracks  and  thus  possibly  result  in  an  engine  ignition  ex¬ 
plosion  accident.  Therefore,  it  is  necessary  to  carry  out 
complete  structural  analysis  of  the  grains.  This  necessitates 
precise  calculations  of  the  pressure  and  strain  values  caused 
by  this  type  of  load.  At  present,  there  already  exist  many 
engineering  calculation  methods:  for  example,  the  approximate 
solutions  [1,2,6]  obtained  when  the  propellant  is  an  incompres¬ 
sible  material  (Poisson  ratio  v=0.5);  or  the  approximate  solu¬ 
tions  [3,5]  obtained  when  the  casing  deformation  is  disregarded; 
in  another  method,  although  the  above  assumptions  are  not  made 
yet  Schapery's  approximate  inversions  were  used  to  obtain 
visocoelastic  solutions.  Although  the  first  two  methods  can 
obtain  simple  calculation  formulas,  yet  the  calculation  results 
tend  to  be  low.  Finally,  because  one  type  of  method  used 
Schapery's  approximate  inversions  directly  induced  from  the 
gradually  varied  process  to  handle  this  type  of  speed  change  of 
ignition,  naturally,  it  was  not  accurate  enough. 

When  deriving  the  viscoelastic  solutions,  this  paper  does 
not  employ  the  above  assumptions  nor  the  approximate  inversions 


and  thus  the  solutions  obtained  are  more  accurate.  Moreover, 
if  we  assume  that  the  propellant's  Poisson  ratio  v=0.5,  we  can 
then  derive  the  calculation  formula  of  the  first  method  from 
these  results;  if  we  disregard  the  casing  deformation  and  take 
E  ^00  ,  we  can  then  derive  the  calculation  formula  of  the 
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second  method  from  these  results.  Therefore,  the  results  of 
this  paper  are  suitable  for  propellants  with  different  Poisson 
ratios  and  also  considers  the  deformation  of  the  casing.  This 
paper  presents  the  final  results  of  two  loading  situations, 
suddenly- applied  pressure  Pi=P?H(t)  and  gradually-applied 
pressure  P^«P?*  [1-exp (-nt)  ] . 

The  following  assumptions  are  made  in  this  paper:  (a)  the 
propellant  is  an  isotropic  and  homogeneous  linear  viscoelastic 
material;  (b)  the  casing  is  a  thin  walled  elastic  small 
cylinder;  (c)  the  Poisson  ratio  v  of  the  propellant  is  a  con¬ 
stant;  (d)  the  propellant  grains  and  casings  are  relatively 
long  and  they  are  in  a  planar  stress  or  planar  strain  state; 

(e)  the  influence  of  the  shape  of  the  inner  hole  is  considered 
in  terms  of  the  stress  concentration  coefficient;  (f)  the  inner 
hole  boundary  does  not  move,  that  is,  we  do- not  consider  the 
ablation  of  the  inner  hole. 

II.  Elastic  Solutions  of  Grain  Stresses  and  Strains 

The  propellant  is  poured  into  a  casing  and  therefore  it 
involves  the  problem  of  analyzing  the  stresses  and  strains  of 
two  joined  cylinders  sustaining  internal  pressure.  See  Fig.  1. 


Fig.  1  Propellant  grains. 

Key:  (a)  Circular  hole  grain;  (b)  Star-shaped  hole  grain. 


We  will  first  study  the  circular  hole  grains.  The  engine 
is  sufficiently  long  and  thus  the  grain  and  casing  are  both  in 
a  planar  stress  or  planar  strain  state.  The  specific  structure 
of  an  engine  possibly  has  the  following  conditions: 

(1)  The  grain  and  casing  are  both  in  a  planar  strain  state 
that  is,  E  z=  £zc=0; 

(2)  The  grain  is  in  a  common  planar  stress  state  and  the 
casing  is  in  a  planar  strain  state,  that  is  cf Z=-P zc=  0. 

It  is  generally  considered  that  an  engine  with  cohesion  of 
the  two  ends  of  the  grain  and  the  casing  is  close  to  the  first 
type  of  condition;  an  engine  with  non-cohesion  (free)  of  the 
two  ends  of  the  grain  is  close  to  the  second  condition.  We 
separately  study  the  solutions  under  both  conditions. 

The  boundary  conditions  of  th  e  two  conditions  are  shown 
in  Table  1. 


j  tt  &  4) 

r  a  a  |  o,(«)  «  -Pi 

r  m  b  J  u  (  b  )  *  u«  <  b  ) 

r  m  e  j  «„<c  )  *  0 

J  a  O,  J  *,  -  0 


*  a  'S) 

o,<  a  )  a  -  p, 

»  <  6 )  *u, (6) 
o„(c) »  0 
O,  *  -  P|.  *».  *  * 


Table  1 


Key:  ( 1)  —  ( 2)  Condition. 


Because  the  deformation  of  the  grain  (inner  cylinder)  is 

constrained  by  the  casing  (outer  cylinder) ,  contact  pressure  P' 

will  be  produced  on  their  connecting  surfaces.  The  contact 

pressure  can  be  determined  based  on  boundary  condition 

u(b)=u  (b)  listed  in  the  table, 
c 

Under  the  first  type  of  condition,  we  obtain 


_ 2(1  -v)q*Pf  _ 

r^+(1  _2v)6.,+£=f;  -f-Cc*+(l-  2V.)W 


As  regards  the  thin  walled  casing,  h/b  ^  1,  and  the  above  form¬ 
ula  can  be  simplified  to 


P' 


_ 2  (  1  -v)P/ _ 

U+C1-2V)W+  ft*-n  Hr  -fr 


(2) 


Under  the  second  type  of  condition,  we  obtain 


p,^__ _ C2at+(61-a*)v;)P, _ 

r  hx  —  nl  ~  IF  V  3  ) 

C(  1  +v)a*+(  l  —  v)6*3  +  ^,— p-(  1  +v.)-|-CCl  -  2v.)6*  +  c1D 

As  regards  the  thin  walled  casing,  the  above  formula  can  be 
simplified  to 

p, - : - »--kx‘-o»w - ___  (4) 

C(i+v)  +  (i-v)W  +  ft"-j)(i_v;)-|-  -|- 


The  last  term  of  the  denominators  of  formulas  ( 1 )  —  ( 4 )  show 
the  influence  of  the  casing's  rigidity.  Often  E/E  1  and 
this  term  is  a  minute  quantity.  Therefore,  we  can  use  the 
Taylor  series  to  expand  the  contact  pressure  formula  into  a 
polynomial  expression  and  disregard  the  minute  quantities  of 
high  order.  For  example,  after  expansion  of  formula  (2) ,  we 
can  obtain 


2  ( 1  -v) 


1  +( 1  -  2  v)x 


TPr 


_ 2  O  _1 _ —(  1  _ v1) 

Cl  +  (  1  — 2v)X‘)1  l+v'-1  h 


E. 


■P, 


=  a,P, -f -P^a^-P" 


(5) 


In  the  formula 


a  -  2  (  1  —  v) 

1  1  +  (  1  -  2  v)xJ 

8 _ ?  (I1-  1  ) 

Pl  C  1  +  (  1  -  2  v)Xl3* 


1  -v 
]  +v 


(1  -V?) 


(6) 

(7) 


In  the  same  way,  after  expanding  formula  (4)  under  the  second 
type  of  condition,  we  can  obtain 


2  +  a* -  l ) v  D  C  2  +  (X»-i  KKl’-.H,  1  _  v?') - 
p  ^T+V+TT^OX5  ‘  ti  +  v+(i  -v)x*ol  u 


E  n 

~E7P/ 


‘  Oj/3  /  Pt 


-P^a.P.-P 


In  the  formula 


-Jr* 


2+(X*-  1  )v 


l+v+(1-v)X>  (10) 

O  C2-KX,-l)v3(XI-l)  n  _ viv 

Ps“  Cl+v  +  (l  -v)x*3l  (  *'  J 

It  can  be  seen  that  under  the  two  conditions,  the  relation¬ 
ship  betwen  contact  pressure  P'  and  internal  pressure  Pi  can  be 
expressed  as 

P>  =  aP,-P'  (11) 

-J7p,  (ii)' 


We  can  see  from  formulas  (7)  and  (10)  that  aside  from  bein'- 
related  to  the  stress-strain  condition,  the  values  of  ei-  and  /S 
are  mainly  determined  by  Poisson  ratios  v,vc  and  /U.  Based  on 
the  assumptions  discussed  in  the  introduction,  Ci.  and  ^  are 
constants  and  when  v=0.5,  <2-  =1, 

After  contact  pressure  P’  is  known,  by  substituting  P*  into 
the  Lame  formula,  we  can  obtain  the  stress  value: 


1/ 1  ol 

\a.  i  - 

p.+  x‘  ( 

1 "  r* . 

)a+  X*  1 

l)  J 

",+  xl-  in 

V  r'  y 

)°— x*“ 

(S-0: 

\p+  x‘  ( 
r,+  x1—  i ' 

,1+—  > 

As  regards  incompressible  material  (v-0.5),  because  &-ml,  the 
coefficients  of  the  first  terms  in  formulas  (12)  and  (13)  are 
equal  to  -1  and  can  thus  be  simplified  to 
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By  substituting  the  known  stress  into  the  strain  formula, 
we  can  then  separately  obtain  the  strain  values  and  the  two 
types  of  conditions. 

Under  the  first  type  of  condition 


Under  the  second  type  of  condition 

e,=  -^r-(a(  — vc.  +  vP,) 

-  1  1  -“«-<* + »>o  -  a)  -r 

a,-vo,+  v/}/) 

=  {5^(0  -  v)(i  -«Xl)  +  (l  +  v)  <1- «>-£-]  +  vj  fi- 

+  xi^~i((1  “v)  +  (l +V)_S"]JT*  (17) 
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we  can  see  from  stress  formulas  (12)  and  (13)  and  strain 
formulas  (14) -(17)  that:  under  the  two  above  mentioned  types 
of  conditions,  c  r  and  £r“  Moreover,  their 

maximum  values  occur  on  the  inner  surface  (r=^)  ,  that  is, 

^max'  (cU  and  £  max'  *6  <*->• 

Below,  we  will  study  the  elastic  solutions  of  the  stresses 
and  strains  of  grains  with  star-shaped  holes. 

When  we  carry  out  a  complete  analysis  of  the  grain  struct¬ 
ures,  we  pay  close  attention  to  the  maximum  values  of  the 
stresses  and  strains.  We  already  pointed  out  above  that  the 
maximum  stresses  and  strains  of  grains  with  circular  holes  occur 
on  the  inner  surface  while  in  star-shaped  grains  the  inner  hole 
star-shaped  angles  have  the  most  serious  stress  concentration. 
Therefore,  the  stresses  and  strains  of  grains  with  star-shaped 
holes  are  maximum  in  the  inner  hole  star-shaped  angle  areas. 

We  record  the  stress's  maximum  value  as  c*(^,t)  and  the 

0 

strain's  maximum  value  as  £*(^,t).  The  upper  symbol  *  re- 
presents  the  value  of  the  grains  with  star-shaped  holes  and  no 
upper  symbol  represents  the  value  of  equivalent  grains  with 
circular  holes.  The  stress  concentration  coefficient  is  de¬ 
fined  as 


k  = 


o?-o; 


0,-0, 

As  regards  the  inner  hole 

o?(  a  )  -o;(  a  )  =k,C<yr(  a  )  -  o,(  a  )) 


(18) 


In  the  formula,  &  *  (  4')=-?^  and  the  stress  difference  of  equiva¬ 
lent  grains  with  circular  holes  is 


or(  a  ) -a,  (  a  )  )=  1  -  a  )P,+P^ 

By  substituting  in  formula  (18),  we  obtain 

o?(«)--p^rr*/C(i -*)P,+P')-P,  (20) 

Below,  we  will  present  maximum  strain  £^(A): 


(19) 
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Under  the  first  type  of  condition 


£«(£,)=  -  1-jr—  [a»(  a  )  — a  )] 


1  -vl 


1  - «)- Jr^']p‘+ «■> 


Under  the  second  type  of  condition 
£»(  a  )  =  — g~(°e(  a  )  —  vo?(  a  )  +vP;] 

— r[-j^r-‘'(1-“)-(1-2v>]p'+-r  ~w^rk'p‘  (22> 

III.  Viscoelastic  Solutions  For  the  Stresses  ana  Strains  of 
Grains 

After  attaining  the  elastic  solution,  by  employing  the  cor¬ 
respondence  principle  in  viscoelastic  mechanics,  we  can  obtain 
the  Laplace  transform  of  the  viscoelastic  solutions  and  then 
obtain  the  viscoelastic  solutions  by  inverse  Laplace  trans¬ 
formation. 

WE  used  the  finite  element  generalized  Maxwell  model  to 
represent  the  mechanical  properties  of  the  propellant  [8]  as 
shown  in  Fig.  2. 


Fig.  2  Generalized  Maxwell  model. 

The  Laplace  transformation  of  the  operator  modulus  is 

E/S _ 


.q(  ll=E  +  y  — i 
?(0“£<+ 


(23) 


9 


The  relaxation  modulus  is  expressed  by  the  Prony  series 


,v 

E(  t  )  =  E.Jr  2  E,  e\p(- 

»  =  1  '  '  ' 


if 


E,=  E(0)=E.+  J  £■/ 


(24) 


■  «  l 


In  the  formula,  N  is  the  Maxwell  element  number;  HZ.  is  the 

.  «  '  ^ 
relaxation  time  f= 


K)' 


We  first  find  the  viscoelastic  solution  P’(t)  of  contact 
pressure  P' .  The  Laplace  transform  of  formula  (11)  is 


?'  ( 0  =  aP.( s )-  p'(  « ) 


(25) 


If  we  use  the  differential  operator  to  express  the  con¬ 
stitutive  equations,  based  on  the  correspondence  principle, 
from  formula  (11) '  we  obtain 


p'(*)=P 


1  Q(  *  ) 

E.  p(s) 


Pi(  *  ) 


(26) 


When  using  the  integral  operator  to  express  the  constitutive 
equations,  then 


P»(s)=P-f-  -gr-s£(s)P-(s) 


(27) 


Below,  we  will  discuss  the  solutions  of  classical  pressur¬ 
ization  equations: 

As  regards  suddenly-applied  pressure  wherein  the  internal 
pressure  is  a  step  function  P^=P?H(t),  at  this  time  P^(s)=P*?/s. 
By  substituting  in  formula  (27),  we  obtain 

P" (s) = /$  —  h—  E  (s)P°,  and  then  seek  the  inverse  and  obtain 

'  n  w .  i 


P'(  o=p  -f- 


(28) 


Therefore 


p'(0=«p;w(o-p-|-  —£pPi 


(29) 
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In  the  formula,  E(t)  can  directly  use  the  measured  value. 

As  regards  the  gradually-applied  pressure,  P . =P? [ 1-exp (-nt) ] , 
and  at  this  time 


^S>=^-T(7-+T,T 


(30) 


We  substitute  formulas  (30)  and  (23)  into  formula  (26)  and 
obtain 


.V  77 

,  n  _  Y  _ V’Pl 

(S  +  "^"r  . -  i  (s  +  «)(*  +-?-)■ 


Seeking  the  inverse,  we  obtain 

^co-P-jp  -fr^s/.o) 


(31) 


In  the  formula 


/■(*)“  i  2  ~trtf—  t  -fr( 

i  »  1 «  v 


*  T'  -e*" 


(32) 


Therefore 


P'  (  t  )  -e-‘)  -  0  -b-  />?/,(  t  )  (33) 

If  we  substitute  formula  (30)  into  formula  (27)  in  the  same  way 
we  can  obtain  formulas  (31)  and  (33).  At  this  time,  auxiliary 
function  fE(t)  is  expressed  as  an  integral  form  and  we  can 
calculate  using  the  measured  E(t)  value: 

£(T)e"VT  (32)' 

By  substituting  derived  pressure  P" (t)  into  stess  formulas 
(12)  and  (13)  we  can  find,  the  viscoelastic  solution  of  the 
stress : 


As  regards  the  suddenly-applied  pressure 


11 


°,(  r  , 

«)_  -A1  f 

:o--£ 

PI 

A:  —  1  1 

.  ■*  J 

E  (  t  ) 

°i(  r  , 

<)_  —  A:  1 

[(■♦* 

P'i 

A1—  1  1 

•P-r 

E(  t  ) 
E. 

1 


)] 

)] 


(- 

(‘+ 


(34) 


(35) 


Therefore 


-  °.(fl>  <) 
Pf  ~  P*, 


CX*(2  a  -  1  )-  n  P  ~ 


E(  < ) 

£. 


(36) 


As  regards  the  gradually-applied  pressure 


+  x^T TT«'> 

a4iri-TE?r((>  +-f-)-tKS-  >))<•  -'*■> 
+T?M1+-£> 


(37) 


(38) 


Therefore 

=  ~  ~Vf  ~ '  =  ’ V~-' 1  CX‘(  2  g~  1 )  ~  1  K  1  -O 

«»> 

Below,  we  will  study  the  viscoelastic  solutions  of  the 
strains.  For  convenience  of  calculation,  we  used  the  finite 
element  generalized  Voigt  model  to  represent  the  mechanical 
properties  of  the  propellant  [8]  as  shown  in  Fig.  3. 
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Fig.  3  Generalized  Voigt  model. 

The  Laplace  transform  of  the  operator  compliance  is 


p(Q 

0(0 


1 


F*+  E  Fl  /  ,  T" \ 
«-»  V+  W 


(40) 


The  creep  compliance  is  expressed  as  the  Prony  series 


FO)  =  F,+  E  Fi  [  1  —exp  (  -j-  )] 


N 


F.-F.+  2  Fi 

I  “  l 


(41) 


In  the  formula,  N  is  the  Voigt  element  number;  .  is  the  delay 

*>i  1 

time  (==— )  . 

Ei 

Under  the  first  type  of  condition,  we  substituted  formula 
(11) ’  into  formulas  (14)  and  (15)  and  used  the  corresponding 
principle  to  obtain 
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-(1 

-  2  v)  (  l  +  v)  /  6* 

1  +  (  1  -  2  v)X! 

as 

7*" 

P.(  *  ) 

(l  * 

2  v)(  1  +v) 

(A1  _ 

l  +  (  1  -  2  v)\* 

V 

PlC  s  ) 

(42) 


If  we  use  the  integral  operator  to  express  the  constitutive 
equations,  then 


wo-  [(1  -2v) 

-t)8TTM,)  <42)< 

w o-  (?"  1 )!?( ’  )W  • >+  -  x^T4  (< 1  '  2  4) 

+4-]»-f  TT'^O  (43)\ 

Pi-P1H (  /  ), 

As  regards  the  suddenly-applied  pressure,  we  substitute  ^ 
P^(s)=P?/s  into  formulas  (42)'  and  (43)',  seek  the  inverse  and 
obtain 

1  )'■ = ■ 1  ~ 2  “> 
-£-]»-r-zr  (4J) 

i?-  - AVc  i -W- 1 )« 1  >+^rH( 1  -  2’} 

<45) 

e„„  -  e«(  a  .  O  _  (  1  -  2  V) (  1  +  *)(X*-  1  )  F(t\ 

PJ  “  P?  1  +  (1  -  2  v)X* 

+i^-o -’’)P-r  E7  (46) 
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As  regards  the  gradually-applied  pressure,  we  substitute 

P.  =P?  [1-exp  (-nt)  ]  and  P.  (s)=P?  -,n  ■  >  into  formulas  (42)  and 

1  1  c  ii  s(stn) 

(43) ,  consider  formula  (40) ,  then  seek  the  inverse  and  obtain 


*Xr,  j  )  (1-2<)(1  +  *)  (  b*  ,  .  'N  _1_  /  (n 

P)  ~  1  +  ( 1  -  2  v;xl  A  r»  )  Et  M  } 

AArM  <  1  -  » ' *>--£-]  »  A  -h 1  (17) 

(v~  -  o  k'’<  <  >+ W[< » —  > 
+-£-)» -f  <«> 

e— -ii (£.  o  -  (1  -  2  y)(  1  +V)  .  1  .  .  v 

PJ  “  PJ"  ■  l+(l-2v)x*'u  m  ' 

,  2  A*(  1  —  v*)  fl  6  1 

+  _____  P  __  )  (49) 


In  the  formula,  auxiliary  function 


If  we  substitute  the  Laplace  transform  of  the  internal  pressure 
into  formulas  (4  2) '  and  (43)  ’  ,  in  the  same  way  we  can  obtain 
formulas  (47),  (48)  and  (49),  only  the  auxiliary  function  is 


/,(  *  )  =  »£>"'  J'  F(i)e”'dT 


A  this  time,  we  can  directly  use  the  measured  value  of  creep 
compliance  (F(t)  for  calculations. 

In  the  same  way,  we  can  find  viscoelastic  solutions  of  the 
strains  under  the  second  type  of  condition.  As  regards  suddenly 
applied  pressure  P^=P°H(t) 


-v)(]  -oX*)-(l  +v)(l-«)-7r-J  +  viF(i) 

+  x^T (<  * -'>-< 1 +’)-?-]  ®T-T7  (51) 

«.('••  O.l'fa  -v)(l  -ain  +  M  +v)(l  -«)  ~7I-]+v}f(  O 
(.  X  1  L 

+  xnrr[(1 -v)  +  (1 +v)^r]P_A"  (52) 

. 1  - 1  -v)(  l  -aX*)  +  (  1  +v)0-  a)X:]+v|  F(  t  ) 


2  X1  q  6  J_ 


(53) 


As  regards  gradually-applied  pressure  P^=P° [1-exp (nt) ] : 


={*irr[(  1  -v)(  1  1  +'’)(>-  a)-pr]+,}'£7/,(  1 > 

5^1- [( i  -v)C  1  -<■!■)+(  l  +■')(!- «)Tr]+viTT/ 

+  i^T[a-v)+0+v,4.]B-|-i7(.-0 

_e«(a»  ii»]  ^  f(l  -v)(l  -aXl)  +  (l  +v)(l—  «)X1]  +  V}£1f  M  ‘  ) 
5  -P  ^  ( x  —  l  (  * 


e,(  r. 

1  f/ 

PI 

“lx1-  1  V 

«.( 

_Q_f  i  f, 

P'i 

“lx1-  1  l 

+  x‘ 

+  X‘-  1 

(54) 


(55) 


2  X1  q  6  (  r  i  _  »-•*> 

——P—  (1  «  ) 


(56) 


In  the  above  formulas,  the  definition  of  auxiliary  function 
f_(t)  is  the  same  as  in  formulas  (50)  and  (50)’.  The  above 

r 

are  viscoelastic  solutions  of  the  stresses  and  strains  of  grains 
with  circular  holes. 

Below,  we  will  study  the  viscoelastic  solutions  of  the 
stresses  and  strains  of  grains  with  star-shaped  holes.  For 
these  types  of  grains,  it  is  only  necessary  to  derive  maximum 
stress  0"*(^,t)  and  maximum  strain  £*(r^,t). 

u  & 

By  substituting  formula  (28)  into  formula  (20)  ,  we  can 


obtain  stress  C*{ctr,  t)  under  suddenly-applied  pressure 


(57) 


By  substituting  formula  (31)  into  formula  (20) ,  we  can 

lly-aj 

b  E. 


obtain  stress  <?*(<2^,  t)  under  gradually-applied  pressure 

G» 


-”S(^  =  _2ii. h,  [(  1  —  a  ) (  1  —  O  +  0  -ft/8 1 C  ° ' 0  1  " 

(53) 

Based  on  elastic  solution  formulas  (21)  and  (22)  for  the 
stress  and  using  the  correspondence  principle,  we  can  find  max¬ 
imum  strain  t)  under  the  two  types  of  conditions: 

& 

Under  the  first  type  of  condition,  as  regards  suddenly- 
applied  pressure  P^=P°H(t) 


(59) 


As  regards  gradually-applied  pressure  P^=P° [l-exp(-nt) ] 


«•(  « .  t)  _  , ,  .  r  2  A*  ,  .  1-2V]  l 

— >{ —  ( 1  ~ v )  l ( 1  - a  >  — r^v-J  e;  m « ) 
+(‘  ~S7(  1  ~e“) 


(60) 


Under  the  second  type  of  condition,  as  regards  the  suddenly- 
applied  pressure 

-*1  V;  i  2  v,  ]  F  ( „ + -££-*' » v  tV 

(61) 

As  regards  the  gradually-applied  pressure 
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i 

I 

s 


> 


x--‘l 

+-0rk^Jt 

Naturally,  if  the  propellant  is  an  incompressible  material 
(v=0.5)  and  Q/-1.,  then  formulas  (57) -(6  2)  above  are  simplified 
into  simple  calculation  formulas  [2]. 

Below,  we  will  use  the  above  mentioned  calculation  formulas 
to  carry  out  calculation  and  analysis  of  specific  propellant 
grains. 

Grain  dimensions:  outer  radius  b=30cm,  inner  radius  a=15.8cm, 
the  star-shaped  angle  number  is  5,  the  star-shaped  angle  coef¬ 
ficient  is  0.8  and  transition  circular  arc  radius  £>=1.5cm,  the 

6  2 

thin  walled  steel  casing:  E  =2x10  kg/cm  ,  Poisson  ratio  v  =1/3, 

c  c 

wall  thickness  h=0.3cm;  the  propellant's  relaxation  modulus 

E (t) =14 . 0+7 . 25exp (-t/0 . 002) +404 , 70exp (-t/0 . 02) +109 . 91exp (-t/0 .2) 

2  2  - 

+500.82exp(-t/2) ,E  =1037kg/cm  ,E  =14kg/cm  ,  Poisson  ratio  v=0.48; 

9  e  o  2 

the  combustion  chamber's  equilibrium  pressure  P.=70kg/cm  and 

transition  pressure  P^=70 [1-exp (-60t) ] . 

1.  Calculation  results  of  suddenly-applied  pressure 
Pi=P°H ( t) =70kg/cm2 : 

We  used  formula  (57)  to  calculate  maximum  stress 
under  two  types  of  conditions  as  shown  in  Table  2: 


f  <  s) 

0 

0.01 

0.03 

0.-05 

007 

0.09 

1  0.10 

1 

og( a  .  1 ) 

tsa® 

0.168 

0.088 

0.085 

-0.025 

-0.041  ' 

-0.050 

-0.053 

-0.165 

1 

«!  a® 

-  0.523 

-  0.565 

-0.606 

1  -  0  623 

-  0.632 

-0-637 

-  0.638 

-0.696 

Table  2 

Key:  { 1 )  —  ( 2 )  Condition. 
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In  order  to  study  the  influence  of  Poisson  ratio  v,  we  also 
made  calculations  of  v=0.5  and  v=0.45.  The  results  are  shown 
in  Fig.  4. 


Fig.  4  t)  under  suddenly-applied  pressure. 

We  used  formulas  (59)  and  (61)  for  the  calculation  of  maximum 
strain  £*(^,t)  under  two  types  of  conditions.  The  results  are 
given  in  Table  3  and  Fig.  5.  When  calculating,  we  used 
F  ( t)  ^  1/E  ( t)  . 


MS)  0  1  001 

0.03 

0.05  1  0.07 

0.09  j  0.10 

1 

te  (  o  .  i  ) 

1»R®|o.»l  x  1 0“ *jo. 90  X  10** 

1.03  x  10-Sjl.ll  x  10'sil.l6  x  10'*  1.19  x  10'*|1. 20X10* 

ll  .83  x  10'* 

P'S 

19K@io.42  x  10'*jo.46  x  10-si 

0.51  x  10"*'0.54  X  10‘S  0.56  x  10‘*|o. 57  x  10**!0. 56  x  10‘S 

0.83  x 10'* 

Table  3 


Key:  ( 1 )  —  ( 2 )  Condition. 


_ \X  y  =  0. 48 

'  v  =  0 . 30' 


0.01  0.03  0.03  0.07  0.03 
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Fig.  5  £g(^,t)  under  suddenly-applied  pressure. 

We  can  see  from  Figs.  4  and  5  that:  under  the  first  type  of 
condition  both  the  stresses  and  strains  are  higher  than  in  the 
second  type  of  situation,  the  stress  decreases  with  the  time, 
the  strain  enlarges  with  the  time  and  they  both  tend  to  be 
limiting  values;  the  influence  of  Poisson  ratio  v  of  the  pro¬ 
pellant  is  very  great  on  the  stress  and  strain  values. 

2.  Calculation  results  of  gradually-applied  pressure 
Pi=70 [l-exp(-60t) ] : 

We  used  formula  (58)  to  calculate  the  maximum  stress  under 
two  types  of  conditions  first  using  formula  (32)  to  calculate 
auxiliary  function  f£(t).  The  calculation  results  are  listed 
in  Table  4  and  Fig.  6. 


0.01 

30.04 


0.05  | 

0.07  | 

Q.r.O 

Oil 

1 

IS. [9  | 

*3-32 

11.40 

to  49 

!  :  2 . 9 1 

0.025  |  - 

0.021 

[  -0.0 to 

1  -0.019 

!  -0.164 

-  0.579  .  - 

0.614 

-0.629 

|  -  0.635 

J  -0.696 

+  V 


/ 


/ 


- 1; 


\ 

PJP « 


• - KihiilE 

— a*!?!  ffi 


0.05  0.07  0.09 


Cl) 

(2) 


Fig.  6  f7*(^/t)  and  P,  (?)  under  gradually  applied  pressure. 

Key:  (1)  Suddenly-applied  pressure;  (2)  Gradually- 
applied  pressure. 

We  used  formula  (50)  to  calculate  auxiliary  function  fp(t) 
and  then  used  formulas  (60)  and  (62)  to  calculate  the  maximum 
strain  under  the  two  types  of  conditions.  The  calculation 
results  are  listed  in  Table  5  and  Fig.  7. 


f<s)  0  1  001 

0.03  0.05  |  0.07  0.09  j  0.11 

1 

/*(<)  |  0  0.507 

1.111  j  1.427  :  1.617  j  1.752  1.859 

3.235 

«o<  a  .  t ) 

WiR®1  0  jo. 39  x  10~3 

0.80X  10'S.0.99x  10"s  1.08  x  10-sjl.l5  x  10  3jl.20x  10‘!: 

1.82  x 10-» 

P° 

MR®'  o  0. 20  x  ]0" 3 

■0.40  x  10‘3j0. 49  x  10"3  0.53  x  10'3  0.56  x  10'3jo.58x  10‘3 

jo. 82  x 10'* 

Table  5 

Key:  ( 1 )  —  ( 2 )  Condition. 
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Fig.  7  £*(<£,,  t)  under  gradually-aoplied  pressure. 

Key:  (1)  Suddenly-applied  pressure;  (2)  Gradually- 
applied  pressure. 

We  can  see  from  Figs.  6  and  7  that:  under  the  first  type  of 
condition,  the  stress  gradually  decreases  from  tensile  stress 
to  compressive  stress  and  under  the  second  type  of  condition 
it  is  compressive  stress  from  beginning  to  end;  under  the 
second  type  of  condition  the  strain  is  markedly  lower  than  in 
the  first  type  of  condition;  the  stresses  and  strains  gradually 
approach  the  value  under  suddenly-applied  pressure  following 
the  time.  When  the  internal  pressure  approaches  the  equili¬ 
brium  pressure,  the  absolute  value  of  the  stress  and  the  value 
of  the  strain  are  both  slightly  lower  than  the  value  of  this 
time  under  suddenly-applied  pressure.  In  engineering  approxi¬ 
mate  calculations,  estimations  of  the  stresses  and  strains  are 
allowable  based  on  calculations  under  suddenly-applied  pressure. 

IV.  Conclusion 

1.  We  can  see  from  the  above  examples  that  some  of  the 
existing  engineering  calculation  methods  assume  that  the  pro¬ 
pellant  is  an  incompressible  material  (v=0.5)  and  the  obtained 
stress  and  strain  values  tend  to  be  low.  In  other  calculation 
methods,  they  easily  disregard  the  terms  with  F'Ec.  We  can  see 
from  formulas  (57) -(62)  that  the  calculation  values  also  tend 
to  be  low,  especially  when  the  casing  material's  elastic  modulus 
Ec  is  relatively  low  (e.g.  aluminum  alloy  and  glass  fiber  rein¬ 
forced  plastic  casings).  However,  use  of  the  method  introduced 
in  this  paper  can  be  used  to  simply,  easily  and  quite  accurately 
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calculate  the  stress  and  strain  values  under  internal  pressure 
loads. 

2.  The  difference  between  the  stresses  and  strains  are  very 
large  under  the  two  types  of  conditions.  Under  the  second  type 
of  condition  (the  casing  is  in  a  planar  strain  state  and  the 
grains  are  in  a  common  planar  stress  state) ,  the  stress  and 
strain  values  are  lowest.  Therefore,  during  calculations,  it  is 
necessary  to  accurately  calculate  which  type  of  stress  state  the 
grains  and  casing  of  the  engine  belong  to.  This  conclusion  also 
shows  us  that  when  designing  an  engine,  we  should  as  much  as 
possible  cause  them  to  be  in  the  second  type  of  condition. 

3.  In  order  to  accurately  calculate  the  stress  and  strain 
values  under  internal  pressure  loading,  it  is  necessary  to  have 
short  time  relaxation  modulus  data  and  accurate  Poisson  ratio 
data.  Both  the  stresses  and  strains  are  very  sensitive  to  the 
Poisson  ratio:  the  smaller  the  v  value,  the  higher  the  stress  and 
strain  values;  when  the  v  value  increases  to  0.5,  the  stress  and 
strain  values  are  lowest.  Moreover,  the  strain  does  not  change 
with  the  time  under  suddenly-applied  pressure. 

After  obtaining  the  stress  and  strain  values  under  internal 
pressure  loading,  the  stress  and  strain  values  are  coordinated 
during  the  solid  cooling,  storage,  environmental  change  and 
transport  etc.  processes  and  we  can  then  further  carry  out  com¬ 
plete  structural  analysis  of  the  propellant  grains. 
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A  STUDY  OF  THE  LOW  FREQUENCY  INSTABILITY  OF  THE  SOLID  PROPELLANT 
ROCKET  ENGINE 

Zhao  Bohua 

Abstract 

This  paper  presents  a  new  concept  of  the  low  frequency  * 
instability  on  the  basis  of  experiments  and  analysis  of  the 
flow  field  of  combustion  gases  in  the  engine.  It  shows  the 
mutual  relationship  between  the  tubular  grain  design  parameters 
(ratio  of  inside  port  to  outside  port  parameter  A/  and  the  com¬ 
bustion  chamber  pressure  p)  and  its  effect  on  low  frequency 
instability.  A  relational  formula  is  derived  through  study  for 
analyzing  the  low  frequency  instability  of  the  tubular  grain 
engine.  The  theoretical  analysis  agrees  with  the  test  results. 
This  paper  has  reference  value  for  selection  of  tubular  grain 
design  parameters  (p  and  7U  )  and  their  inhibition  of  low  fre¬ 
quency  instability. 

I.  Introduction 

Low  frequency  instability  such  as  low  frequency  oscillation 
and  surge  is  generally  a  nonacoustic  oscillation  problem.  The 
p-t  test  curve  of  the  pressure  of  commonly  seen  low  frequency 
oscillation  combustion  and  surge  which  changes  with  time  is 
shown  in  Fig.  1. 


CO  w 


Fig.  1  p-t  curve  of  low  frequency  oscillation  and  surge. 

In  the  figure,  (a)  is  the.  increased  amplitude  low  frequency 
oscillation  in  the  engine  due  to  certain  factors  having  very 
small  vibration.  At  this  time,  amplitude  lengthened  index 
s  >0.  (b)  is  the  decreased  amplitude  low  frequency  oscilla¬ 

tion  combustion  in  the  engine  due  to  certain  factors  producing 
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certain  amplitude  fluctuations.  At  this  time,  amplitude  atten¬ 
uation  index  s  K.  0  which  shows  that  the  damping  in  the  engine 
causes  the  oscillation  to  gradually  weaken  and  finally  tend 
towards  stable  combustion.  (c)  is  the  p-t  curve  of  equal 
amplitude  low  frequency  oscillation  combustion  when  s=0. 

(d)  is  the  surge's  p-t  curve.  There  are  two  types  of  low  fre¬ 
quency  instability,  linear  and  nonlinear.  This  paper  studies 
the  problem  of  linear  low  frequency  instability  with  oscilla¬ 
tion  frequency  of  about  100Hz. 

As  regards  the  tubular  grain  engine,  when  the  operating  time 
of  the  engine  is  short  and  there  are  strict  requirements, 

X  .  . 

A*  ^  l  is  often  used  and  at  this  time  relatively  large 

outside 

loading  density  can  be  obtained.  Tests  show  that  if  the  operat¬ 
ing  pressure  of  the  combustion  chamber  is  greater  than  the 
critical  pressure  value  of  this  propellant  and  the  A*  value  of 
the  tubular  grain  is  smaller  than  a  certain  A*,  at  the  low  of 
-40°C,  low  frequency  instability  combustion  will  occur. 

This  type  of  low  frequency  instability  will  cause  the  thrust 
to  have  relatively  large  fluctuations  and  changes.  It  also 
often  causes  rockets  and  missiles  to  have  poor  firing  concentra¬ 
tion  at  low  temperatures  and  even  the  dropping  of  bombs.  When 
the  low  frequency  oscillation  frequency  is  equal  to  the  natural 
oscillation  frequency  of  the  rocket  or  missile  system,  this  can 
result  in  the  mechanical  oscillation  of  the  entire  missile 
system  or  components  causing  the  missile  to  be  unable  to  operate 
normally.  In  certain  strategic  missile  engines  as  well  as  in 
aerospace  applications,  because  it  is  necessary  to  use  the  ad¬ 
vantages  of  low  combustion  chamber  pressure  and  to  raise  the 
specific  impulse  of  the  weight  of  the  engine,  study  and  resolu¬ 
tion  of  this  problem  has  important  and  real  significance. 

II.  The  Problem  of  Low  Frequency  *  Instability 

If  there  is  no  heat  loss  in  the  solid  phase  of  the  propell¬ 
ant,  its  heat  conduction  control  equation  is 
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a  = 


In  the  formula, 


p,c,’ 


separately  indicate  the 


density  specific  heat  and  heat  conduction  coefficients  of  the 
solid  propellant  and  r  is  the  combustion  rate  of  the  propel¬ 
lant.  The  source  point  of  the  x  axis  is  taken  on  the  surface 
of  the  propellant  as  shown  in  Fig.  2. 


Fig.  2  Chart  of  the  temperature  distribution. 

The  boundary  conditions  are 

X=—  ,  T=T^ 

X=-0 ,  T=T 

s 

P  T  =H_  (temperature  gradient  when  x=-0) 

ax  s 

Based  on  the  above  boundary  conditions,  by  solving  partial 
differential  equation  (1) ,  we  obtain  the  temperature  distribu¬ 
tion  formula  of  the  solid  propellant  grains  when  in  a  steady 
state 

7=(^)”'exp[{~^)x)+T'  (2) 

When  x=-0  and  T=T  are  substituted  into  formula  (2) ,  we 

s 

obtain  the  grain  surface  temperature  formula 


When  in  a  steady  state,  the  energy  equilibrium  relation 
in  the  solid  propellant's  grain  surface  area  can  be  written  as 

KH,  =  KtH,.+  P,rQ  (4) 

In  the  formula,  Q  is  the  surface  phase  change  heat. 

When  formula  (4)  is  multiplied  by  r/KHg  and  substituted  into 
formula  (3),  we  obtain  the  expression  of  the  grain's  steady 
state  combustion  rate. 

r=5  if(T  _ T)  C  5  ) 

Above,  we  obtained  temperature  T.  of  the  solid  propellant 
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surface  under  a  steady  state  condition  as  well  as  the  theoret¬ 
ical  expression  of  the  grain's  combustion  rate  r.  If  there  are 
certain  factors  which  cause  the  grain's  surface  temperature  to 
change,  it  then  necessarily  directly  influences  the  change  of 
the  grain's  combustion  rate  and  the  change  of  the  combustion 
rate  will  cause  changes  in  the  combution  chamber  pressure. 

As  regards  the  tubular  grain  engine  with  A  K.  1,  because 
flow  rates  of  the  combustion  gases  of  the  inside  port  and 
outside  port  are  not  equal,  the  parallel  gas  flow  of  the 
outside  port  must  be  divided  into  two  gas  flows  by  the  stagna¬ 
tion  point  and  the  position  of  the  stagnation  point  is  related 
to  the  size  etc.  of  the  A  value.  Figure  3  shows  the  model  of 
the  combustion  gas  flow  field  in  this  type  of  engine.  Because 
the  combustion  gas  flows  of  the  outer  surface  of  the  solid 
propellant's  tuburlar  grain  flows  backwards,  under  certain 
conditions,  it  possibly  causes  disturbance  of  the  combustion 
gas  flow  in  the  vicinity  of  the  stagnation  point  and  the  forma¬ 
tion  of  helicon  waves.  It  is  thus  necessary  to  irritate  the 
grain  surface  in  varying  degrees.  When  the  ratio  of  inside 
port  to  outside  port  parameter  A  of  a  single  tubular  grain 
engine  with  double-base  propellant  is  0.6,  at  -40°C,  the 


combustion  chamber's  equilibrium  pressure  is  7.85MPa  and  the 
test  curve  of  the  low  frequency  oscillation  pressure  which 
changes  with  time  is  as  shown  in  Fig.  4.  When  observed  from 
the  finally  combusted  grain,  there  are  grain  spots  and  ripples 
as  well  as  the  formation  of  inverted  cones  on  the  grain's 
outer  surface  of  the  grain's  front  end  (far  from  the  nozzle 
end) .  The  final  test  proved  the  accuracy  of  the  hypothesis  of 
the  combustion  gas  flow  field  model  in  the  engine. 

V7e  can  therefore  simply  describe  this  type  of  low  frequency 
oscillation  combustion  as  follows:  under  certain  pressure  condi¬ 
tions,  when  the  7^  value  is  smaller  than  a  certain  7V  *  value, 
because  of  the  backward  flow  of  the  combustion  gas  flow,  a  very 
small  disturbance  of  the  combustion  gas  flow  will  be  caused  near 
the  stagnation  point  of  the  outer  surface  of  the  grain.  The 
disturbance  or  oscillation  of  the  gas  phase  in  the  engine  can 
cause  periodic  changes  in  the  heat  conduction  and  temperature 
and  the  decomposition  speed  of  the  solid  phase  can  be  taken  as 
the  function  of  the  periodic  changes  of  the  grain  's  surface 
temperature.  Based  on  the  feedback  mechanism,  the  fluctuation 
of  the  combustion  rate  uses  the  fluctuation’ of  the  pressure 
caused  by  the  combustion  chamber's  transfer  function  and  the 
pressure  fluctuation  uses  the  feedback  of  the  combustion's  trans¬ 
fer  function  and  the  resulting  fluctuation  of  the  combustion  rate. 
As  a  result,  low  frequency  oscillation  combustion  can  be  produced 
in  the  solid  propellant  rocket  engine. 


Fig.  3  Schematic  of  combustion  gas  flow  field  in  engine  when 
7^<  l. 

Key:  (1)  Front  end;  (2)  Stagnation  point;  (3)  Rear  end. 


Fig.  4  Curve  of  low  frequency  oscillation  pressure. 

We  know  based  on  the  energy  conservation  equation  that  in 
the  combustion  chamber  the  combustion  gas  produced  from  the 
surface  of  the  propellant  grain  is  equal  to  the  mass  flow 
flowing  passed  the  nozzle  as  well  as  the  sum  of  the  changes  of 
the  combustion  gas  mass  in  the  combustion  chamber.  That  is 


AtPpT- 


CD  . 
'~~g~A,P 


d  (p^.)  _  c 
dt 


Because 


then 


L_  At. 

dt  gRT t  dt 


(pf~ P,)Atr-  g°  A'P  +  gRTt  dt 

In  the  formulas,  indicates  the  combustion  surface  of  the 
grain;  p ^  and  separately  indicate  the  density  and  tempera¬ 
ture  of  the  combustion  gas;  R  is  the  gas  constant;  t  is  the  time; 
C_  is  the  weight  flow  coefficient,  V  is  the  free  volume  of  the 
combustion  chamber;  Afc  is  the  nozzle  throat  area;  g  is  the 
gravitational  acceleration. 

If  low  frequency  oscillation  combustion  is  produced  in  the 
engine  as  shown  in  Fig.  4,  we  assume  the  pressure  amplitude  of 
the  combustion  chamber  assumes  a  sinusoidal  shape.  That  is 

£=p  —  p0e"sir<io/  (  7  ) 
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If  the  phase  angle  of  the  combustion  phase  of  the  pro¬ 
pellant  grain  to  the  pressure  low  frequency  oscillation  lead 
is  ,  then 


r  =<jp’=  a  (_p- p„e,'sin(La>t+  P)3* 

If  we  expand  the  right  end  series  of  the  above  formula  and  dis¬ 
regard  the  higher  term,  then 

r  a  (p*— np„p"‘ V'sin(<of  +  (3  )]  (  8  ) 

In  the  formula,  p  is  the  pressure  when  there  is  steady  state 
combustion;  pQ  is  the  assigned  initial  moment  pressure  ampli¬ 
tude;  a  is  the  combustion  rate  coefficient;  n  is  the  pressure 
index;  Gs  is  the  oscillation  angle  frequency;  s  is  the  amplitude 
lengthening  index  or  attenuation  index. 

When  s  ^  0,  this  shows  that  the  pressure  amplitude  lengthens 
with  time  based  on  the  index  law;  when  s  D,  this  shows  that 
the  pressure  amplitude  attenuates  based  on  the  index  law. .  At 
this  time,  s  indicates  the  size  of  the  low  frequency  oscilla¬ 
tion  damping;  when  s=0 ,  this  indicates  equal  amplitude  low  fre¬ 
quency  oscillation. 

By  differentiating  formula  (7) ,  we  obtain 

—  (p,e,,<»cos<i>/+  p„e"sshKi>0  (9) 

Because  combustion  gas  density  P  and  propellant  density 
P  can  be  disregarded  and  not  calculated  by  substituting  form¬ 
ulas  (7) ,  (8)  and  (9)  into  formula  (6) ,  we  obtain 


•  i 

-?—tr—anfytpm''e''s\n(<i>t+  P)) 

T 

=  p-/yr"siT>a>/ />,*"(  ®  cos  of  +  *  sinco<)5 
— |—  [unp0e'’p"'sin((i>t+  P  )) 

-  p.‘"  1  »''“'+rafc.c““s“'+ssi''“°! 

Letting  * 


’-cfc-  Wl,**u>  <10) 

C=anp„e"p-*  (2) 


Key:  (1)  Time  constant;  (2)  Combustion  rate  disturbance 
function;  (3)  Pressure  disturbance  function. 


Then 


R 

'<D 


sin(co/  +  p  )  a  (  i  -fsT)sinof +  CDtcoscof 


=  (  1  +5T)J]1/SsiTlC®^+  P  ) 


•  < 


In  the  formula 


cos  P  = 


1  +5T 

C(»T)*+(1  +  ST)*}"*  ’ 


sin  p  = 


_ cor 


After  simplification,  we  obtain 

*J1  - _ 1 _ 

t/f  ”  C(COT)l+(l  +JT)0'/1  "  } 

Formula  (11)  indicates  the  relational  formula  of  the  combustion 

chamber's  transfer  function  of  the  pressure  fluctuation  caused 

—2 

by  the  combustion  rate  fluctuation.  By  multiplying  r  /4a  by 
/2Z  in  the  formula,  we  obtain  dimensionless  time  constant 
TZ  -7  ^/4a  and  the  critical  time  constant  if  there  is  no  low 
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frequency  oscillation  combustion  is  (/‘'n)cr'  that  is 


(0.r  = 


oW. 

4aC oRT  gj4t 


This  dimensionless  critical  time  constant  reflects  the  proper¬ 
ties  of  the  solid  propellant,  the  special  features  of  the 
engine's  structure,  the  design  parameters  of  the  grain  and 
other  major  factors  influencing  low  frequency  instability.  It 
is  a  comprehensive  parameter. 


We  let 


V,  Di-P'  +  d^L 

A,  ~  4 A, 

„  xi1}  D\-Dl 


Key:  (1)  Inside;  (2)  Outside. 

In  the  formulas,  is  the  inside  diameter  of  the  combustion 
chamber;  D,  d  and  L  are  separately  the  outside  diameter,  inside 
diameter  and  length  of  the  grain;  ■t  is  the  correction  coeffic¬ 
ient  taking  into  consideration  the  ignition  space  and  grain 
supporting  space  in  relation  to  free  volume  Vc;  xQUtside  is  the 
inside  port  parameter;  /L.  *  is  the  characteristic  ratio  of 
inside  port  to  outside-port. 

We  substitute  formula  (13)  into  formula  (14)  and  letting 
S^=7TDdL  and  f)  =^,  we  then  obtain 

v,  _  +S/(X*+  *0  (15) 

At  4  At 

By  substituting  formula  (15)  into  formula  (12),  we  obtain 


16aC0rr^,(T.),f  ,  » ,  *1"  _  n 


We  can  rewrite  formula  (16)  as  another  equivalent  form,  that  is, 
the  lowest  critical  pressure  p*  which  does  not  produce  low  fre 
quency  oscillation  combustion 
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(16)' 


\ZaCDRTtA,(T*) 

iji5’,a1(X*+  *1  ) 


Formula  (16)  shows  that  the  index  relationship  of  /b  *  and  p* 
have  a  (-2n)  linear  relationship  to  the  given  solid  propellant 
and  tubular  grain  engine.  In  order  to  guarantee  that  low 
frequency  oscillation  combustion  does  not  occur,  the  smaller  the 
7b*,  the  higher  the  required  critical  pressure  value  p* . 


III.  Analysis  of  Test  Results 


When  a  certain  single  double-base  tubular  grain  engine 
at  -40°C,  the  equilibrium  pressure  of  the  combustion 
chamber  is  7.85MPa  and  the  engine  has  stable  combustion.  Under 
other  constant  conditions,  we  carried  out  several  types  of  ex¬ 
periments  with  /V  values  less  than  1  and  when  the  7 b  value  was 
less  than  a  certain  *  value,  then  low  frequency  oscillation 
combustion  occurred.  The  measured  test  curve  is  shown  in 
Fig.  5. 


(1) 


(2)  (3) 


Fig.  5  p-t  test  curve  when  the  double-base  propellant  is  at 

-40°C,  the  pressure  is  7.85MPa  and  there  are  different 
7U  values. 

Key:  ( 1 )  —  ( 5 )  Outside. 
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By  analyzing  the  test  results,  we  know  that: 

1.  As  regards  a  tubular  grain  engine  with  A^l,  when  the 

A  value  is  less  than  a  certain  7i*  value,  at  a  low  temperature 
of  -40 °C,  low  frequency  oscillation  combustion  first  occurs. 
Different  trademark  propellants  have  different  A*  values  of  low 
frequency  oscillation  combustion  under  the  same  combustion 
chamber  pressure  conditions.  The  larger  the  low  frequency 
oscillation  amplitude  becomes,  the  lower  the  average  pressure  and 
when  the  pressure  is  lower  than  a  certain  value,  combustion  is 
extinguished  or  surges. 

2.  The  smaller  the  A  value,  the  greater  the  pressure  decrease 
loss  A p  and  the  faster  the  time  the  grain  cuts  off  combustion. 

The  test  curve  of  double-base  propellant  A  p/p  which  changes  with 

7L  is  shown  in  Fig.  6. 


Fig.  6  Test  curve  of  double-base  propellant  A  p/p  which  changes 
with  A  at  -40°C. 

Key:  (1)  Single  tubular  grain. 

3.  The  measured  low  frequency  oscillation  frequency  is  about 
100Hz,  the  low  frequency  oscillation  frequency  is  almost  unrelated 
to  the  size  of  A*  and  linearly  increases  with  the  base  pressure 
(i.e.  the  lowest  pressure  between  the  adjacent  pressure  peaks). 

The  amplitude  of  low  frequency  oscillation  combustion  generally 
decreases  with  the  increases  of  the  base  pressure  and  thus  the 
amplitude  generally  decreases  with  the  increases  of  the  frequency. 
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4.  The  smaller  the  TV*  value,  the  higher  the  lowest  critical 
pressure  value  which  does  not  produce  low  frequency  oscillation 
combustion.  It  was  observed  in  the  tests  that  after  low  frequency 
oscillation  combustion  occurred,  if  we  suitably  raised  the 
combustion  chamber  pressure  or  increased  the  TV  value,  this  could 
eliminate  the  low  frequency  oscillation  combustion.  The  partial 
test  results  of  a  double-base  propellant  tubular  grain  engine 
wherein  TV*  changes  with  p*  at  -40 °C  are  shown  in  Fig.  7.  The 
test  results  agree  with  theoretical  formula  (16). 


Fig.  7  Test  curve  of  A*  -  p*  of  double  base  propellant  tubular 
grain  engine  at  -40°C. 

Key:  (1)  Low  frequency  oscillation  combustion  area; 

(2)  Stable  combustion  area. 

IV.  Major  Factors  Influencing  Low  Frequency  Instability 
Low  frequency  instability  has  been  observed  when  using 
common  double  base  propellants,  double-base  propellants  containing 
lead  salt  and  metal  additives  as  well  as  in  other  types  of 
propellants;  low  frequency  unstable  combustion  also  occurred  in 
engines  with  different  types  of  grains  such  as  end-burning  grain, 
tubular  grain,  internal-burning  star-shaped  grain  etc.  The 
factors  influencing  low  frequency  instability  are  complex  and  we 
can  know  from  formula  (16)  and  analysis  of  the  test  results  that 
the  major  factors  influencing  solid  propellant  rocket  engine 
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instability  are  the  composition  and  properties  of  the  solid 
propellant,  the  grain's  initial  temperature,  the  structural 
characteristics  of  the  engine  and  the  grain's  design  parameters 
etc.  We  will  present  a  simple  analysis  below. 

1.  Composition  and  Properties  of  the  Solid  Propellants 

The  low  frequency  instability  and  critical  pressure  values 

of  solid  propellants  with  different  compositions  are  different. 
Generally,  double-base  propellants  which  do  not  contain  metallic 
powder  but  have  organic  lead  salt  content  and  uniform  distribu¬ 
tion  have  marked  effects  on  low  frequency  instability  and  the 
lower  the  content  the  worse  its  low  frequency  instability. 

Because  organic  lead  salt  is  a  type  of  combustion  catalyst,  under 
low  pressure  it  has  important  regulating  effects  on  the  combustion 
rate  of  solid  propellants.  The  catalytic  effects  are  also  related 
to  the  size  of  the  organic  salt  so  that  the  smaller  the  size,  the 
more  uniform  the  distribution  in  the  absorbed  grain  and  the 
better  the  catalytic  effects. 

The  addition  of  aluminum  powder  or  another  metal  as  well  as 
metal  oxide  is  a  good  method  for  inhibiting  low  frequency 
oscillation  combustion  yet  it  also  has  its  disadvantages  in  its 
effects  on  low  frequency  instability.  Because  the  reaction  heat 
of  metallic  oxide  is  quite  high,  under  low  temperatures,  it  can 
help  solid  propellant  thermal  decomposition  and  accelerate  the 
combustion  process.  Aside  from  this,  light  metals  are  effective 
for  inhibiting  the  pressure  disturbance  of  the  combustion  process. 
If  the  grain  size  and  contents  of  the  metallic  powder  in  the 
solid  propelllant  ate  unsuitable,  low  frequency  oscillation 
combustion  can  occur  during  combustion  because  of  the  periodic 
falling  off  of  the  accumulated  metal  from  the  propellant  surface. 

2.  Initial  Temperature  of  Grain 

Given  that  other  conditions  are  the  same,  it  is  much  easier 
for  low  frequency  unstable  combustion  to  occur  when  the  grain  is 


at  a  low  temperature  than  when  at  normal  and  high  temperatures. 
This  is  because  when  at  high  and  normal  temperatures,  the 
combustion  chamber  pressure  is  generally  larger  than  the  pressure 
at  low  temperatures.  The  higher  the  pressure  the  larger  the 
cumbustion  gas  density  and  the  larger  the  damping  effect  on  the 
initial  disturbance  in  the  combustion  gases.  Raising  the  pressure 
has  important  effects  on  inhibiting  the  low  frequency  oscillation 
combustion  caused  by  pressure  disturbance.  Therefore,  for  a 
specific  engine,  it  is  only  necessary  that  it  not  have  low 
frequency  oscillation  combustion  at  low  temperatures  and  then  it 
will  have  stable  combustion  at  normal  and  high  temperatures. 

3.  Structural  Properties  of  the  Engine  and  the  Grain's 
Design  Parameters 

The  influence  of  the  structure  of  the  engine  on  low  frequency 
instability  is  very  complex.  For  example,  when  other  conditions 
are  constant,  the  low  frequency  stability  performance  of  reverse 
flow  multiple  nozzle  single  grain  engines  is  better  than  that  of 
direct  flow  single  grain  engines.  This  is  also  to  say  that  under 
tne  same  pressure  conditions,  the  A>  *  of  the  former  is  smaller 
than  that  of  the  latter  or  that  under  the  same  /V*  value  condi¬ 
tions  the  allowable  critical  pressure  p*  value  of  the  former  is 
lower  than  that  of  the  latter.  Because  of  the  complexity  of  the 
problem,  at  present,  there  is  still  no  way  of  theoretically 
calculating  the  influence  of  the  structural  properties. 

Different  grain  types  have  different  influences  on  low 
frequency  combustion  instability.  Tests  show  that  the  critical 
pressure  of  end-burning  grain  engines  is  much  lower  than  that  of 
other  grain  type  engines  because  the  combustion  of  end-burning 
grain  is  not  influenced  by  the  flow  rate  of  combustion  gas  and 
the  geometric  parameters  of  the  grain  shape.  Therefore,  under 
the  same  pressure  conditions,  similar  propellants  are  stable  in 
end-burning  grain  engines.  However,  they  can  become  unstable  in 
tubular  or  star-shaped  grain  engines.  This  is  then  the  influence 
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of  the  grain  type's  geometric  parameters  on  low  frequency  in¬ 
stability.  Internal-burning  star-shaped  grain  also  possibly 
has  rotation  similar  to  a  vortex  flow  in  a  star-shaped  cavity 
as  shown  in  Fig.  8. 


Fig.  8  Vortex  flow  of  star-shaped  cavity. 

Under  certain  conditions,  this  type  of  rotation  must  in  varying 
degrees  puncture  the  inside  surface  of  the  grain  and  cause 
periodic  fluctuations  of  the  solid  propellant  grain's  surface 
temperature  and  combustion  rate.  This  can  bring  about  low 
frequency  unstable  combustion  and  the  star-shaped  grain’s 
geometric  parameters  are  naturally  also  an  important  factor 
causing  low  frequency  instability. 

The  larger  the  grain's  port  parameter  x,  the  higher  the 
critical  pressure  value;  under  certain  pressure  conditions,  the 
larger  the  x  value  the  higher  the  allowable  A*  value,  yet  its 
level  of  influence  is  not  great. 

To  sum  up,  the  critical  pressure  needed  to  guarantee  that  a 
tubular  grain  engine  does  not  have  low  frequency  instability  is 
generally  related  to  the  solid  propellant’s  type  (SP) ,  the 
lowest  operating  temperature  (-T) ,  the  grain's  design  parameters 
and  the  structural  properties  of  the  engine  (CC)  etc.  Therefore, 

/>*-/,( SP,  - T ,  *,  X*.  CC)  (17) 

As  regards  given  solid  propellants,  the  lowest  critical 
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pressure  for  tubular  grain  engines  with  x  x*  and  A  K  1  not 
to  have  low  frequency  unstable  combustion  can  be  simplified 
and  expressed  as 

'  '’-’-'‘IB-)”]  a8) 

Each  type  of  solid  propellant  can  obtain  a  p*T  curve  which 
changes  with  A*  from  tests.  If  low  frequency  unstable  combustion 
occurs  during  the  development  process  of  the  product,  we  can 
know  from  formula  (16)’  or  formula  (18)  that  it  is  only  necessary 
to  suitably  raise  the  value  or  pressure  p  value  to  be  able  to 
completely  inhibit  low  frequency  unstable  combustion. 
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IGNITION  AND  INITIAL  PRESSURE  PEAK  OF  SMALL  SOLID  PROPELLANT 
ROCKET  ENGINES 
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Chen  Hongzhang  and  Jia  Qingying 
Abstract 

This  paper  mainly  studies  the  causes  of  the  initial 
pressure  peak  of  small  solid  propellant  rocket  engines  which 
form  at  high  temperatures,  proposes  a  specific  method  for 
predicting  the  ignition  pressure  peak  and  discusses  the  specific 
measures  for  eliminating  or  reducing  the  initial  pressure  peak 
at  high  temperatures.  These  studies  are  of  practical  signifi¬ 
cance  to  further  improvements  of  the  design  of  solid  rocket 
engines  and  raising  their  overall  performances. 
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correction  coefficient  of  flow  of  engine 
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:  grain's  combustion  surface  area 
:  grain's  initial  combustion  surface  area 
:  propellant's  combustion  rate 
:  ignition  time  of  propellant 
:  time  the  maximum  ignition  pressure  occurs 
:  time  the  grain  combustion  finishes 
:  time  corresponding  to  P  ^ 

:  any  given  time  in  operation  of  combustion  chamber 
:  production  quantity  of  combusiton  gases  of  propellant 
:  pressure  peak  ratio 

:  production  quantity  of  black  powder  gases 
:  production  rate  of  propellant's  combustion  gases 
:  production  rate  of  black  powder  gases 
:  initial  free  volume  of  combustion  chamber 
:  combustion  rate  coefficient  of  propellant 
:  combustion  rate  pressure  index  of  propellant 
:  length  of  combustion  chamber 
:  propellant's  grain  length 
:  inside  diameter  of  combustion  chamber 
:  outside  diameter  of  grain 
:  inside  diameter  of  grain 
:  end  free  port  area  of  grain 

:  ratio  of  initial  combustion  surface  area  and  grain's 
end  free  port 

:  meat  thickness  at  any  given  time  grain  burns  out 

:  meat  thickness  of  burn  out  when  combined  combustion 
period  is  completed 

:  mass  fraction 

:  specific  heat  ratio  of  mixed  gases 
:  specific  heat  ratio  of  black  powder  gases 
:  specific  heat  ratio  of  propellant's  gases 


I.  The  Combustion  Properties  of  Black  Powder  and  the  Ignition 
State  of  the  Engine 

1.  Combustion  Properties  of  Black  Powder 

Based  on  experimental  research  on  the  combustion  properties 


of  black  powder,  we  consider  that: 

(A)  The  combustion  rate  of  black  powder  is  very  fast  and  when 
the  black  powder  flows  passed  the  nozzle  and  there  is  not  a  great 
loss  (aside  from  when  black  powder  is  placed  on  one  end  of  the 
nozzle),  the  maximum  ignition  pressure  point  in  the  engine's 
ignition  pressure  curve  is  the  finishing  point  of  the  combustion 
of  the  black  powder. 

(B)  There  is  very  little  influence  of  environmental  pressure 
on  black  powder  combustion. 

2.  Ignition  State  of  the  Engine 

(A)  Optimum  Ignition  State 

The  optimum  ignition  state  is  when  the  ignition  pressure  in 
the  engine  reaches  maximum  value,  the  initial  combustion  surface 
of  the  entire  propellant  is  ignited.  Naturally,  under  these 
circumstances,  the  pressure  establishment  process  of  the  com¬ 
bustion  chamber  is  ideal  and  it  can  produce  an  ignition  pressure 
peak. 

(B)  Early  Ignition  State 

Early  ignition  of  an  engine  points  to  the  complete  igniting 
of  the  propellant's  combustion  surface  occurring  far  before  the 
ignition  pressure  reaches  the  maximum  ignition  pressure.  In 
this  type  of  situation,  there  exists  in  the  pressure  establish¬ 
ment  process  of  the  combustion  chamber  the  combined  combustion 
period  of  black  powder  and  propellant.  Its  existence  is  mainly 
due  to  the  formation  of  the  ignition  pressure  peak.  Because  the 
propellant  requires  different  ignition  energy  under  different 
initial  temperatures  and  the  ignition  powder  quantity  generally 
used  for  the  engine  is  difficult  to  change,  the  ignition  of  the 
small  solid  engine  under  high  temperatures  is  in  a  serious 
early  ignition  state  which  is  difficult  to  avoid.  Figure  1 
shows  the  pressure  time  curve  of  the  early  ignition  engine. 


Fig.  1  Pressure  time  curve  of  an  engine  in  early  ignition  state. 

II.  Prediction  of  Ignition  Pressure  Peak 

1.  Basic  Formulas  for  Calculating  Pressure  in  Combined 
Combustion  Period 

The  process  of  pressure  change  in  the  combustion  chamber  of  a 
solid  rocket  actually  shows  the  contradictory  process  between 
the  production  and  elimination  of  combustion  gases.  This  paper 
begins  from  actual  engineering  applications,  disregards  some 
factors  of  secondary  importance  and  proposes  the  following 
assumptions  for  the  operating  process  of  the  combined  combustion 
period  of  the  solid  rocket  engine: 

(A)  The  temperature  and  pressure  of  the  combustion  gases  in 
the  entire  combustion  chamber's  free  volume  are  uniform. 

(B)  The  combustion  gases  in  the  combustion  chamber  can  be 
used  for  ideal  gas  processing. 

(C)  The  ignition  of  the  entire  initial  combustion  surface 
occurs  simultaneously. 

(D)  The  exhaust  of  the  nozzle  is  only  carried  out  at  the 
speed  of  sound. 

Based  on  the  principle  of  the  conservation  of  mass  in  the 
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combined  combustion  period,  the  production  quantity  of  combustion 
gas  should  be  equal  to  the  combustion  gas  flow  passed  the  nozzle 
in  addition  to  the  change  rate  of  the  stored  quantity  of  combus¬ 
tion  gases  in  the  combustion  chamber.  We  use  the  following 
equation  to  express  it 


CaF(k)A,P.  .  „  Vt  dP.  P.Vt  d(C,f, r) 

P,su  +  mit=  *^——  +  PkSu+-cj7  d(  {Cjky  dt 


(i) 


After  transposition  and  disregarding 
high  order,  we  obtain 


-iall  quantity  of 


dp.  c.t,  r .  . .  cj?(h)A.p. ,  py,_  t-,  i 

if - /C7,  -+lcjjr  dt  1  '■2> 

Here,  the  focus  of  analysis  is  on  the  expression  of  the 

propellant's  combustion.  It  is  well  known  that  under  non-corrosion 

effects  (e.g.  end  combustion) ,  its  combusiton  rate  formula  can  be 

expressed  by  usu-^p”.  However,  in  an  engine  with  serious  corrosive 

combustion,  the  combustion  rate  on  different  surfaces  of  the 

propellant  is  the  function  of  the  pressure  and  is  also  related  to 

the  gas  flow  rate  etc.  in  this  area.  Tests  show,  however,  that 

as  regards  small  solid  engines,  especially  engines  with  end 

combustion  and  side  combustion  engines  under  L/D  ^  4  ^  A  /A 

P  ^ 

conditions,  we  need  not  consider  the  influence  of  corrosion 
effects  on  the  combustion  rate.  When  carrying  out  combustion 
chamber  P-t  curve  pressure  calculations,  we  can  still  use  the 
previously  given  combustion  rate  formula.  When  carrying  out 
integration  operations  in  a  relatively  large  pressure  range,  we 
should  select  a  corresponding  combustion  rate  formula  with  a 
different  pressure  range  in  a  different  pressure  area  of  the 
integration. 


2.  Calculation  of  the  Combustion  Chamber's  Internal  Pressure 
in  the  Combined  Combustion  Period 


45 


(A)  Determination  of  the  Thermophysical  Properties  of  the 
Mixed  Gas 

in  order  to  approximately  determine  the  main  thermophysical 
properties  of  the  mixed  gas,  we  now  introduce  a  simplified 
formula  to  determine  mixed  gas  powder  power  f ^ [ 1 ] 

f  a-P>//.+P/ 

. 

In  the  formula 


**«(i -3)^+0*# 

6  = — ^ — - 
m,  +  mia 

From  the  above  formula,  we  can  easily  determine  the  values  of 
powder  power  f,  and  specific  heat  ratio  function  F(k)  of  the 

jC 

mixed  gas  at  different  moments  during  the  combined  combustion 
period. 

(B)  Determination  of  the  Gas  Production  Rate  of  Black  Powder 
Figure  2  is  the  ignition  pressure  curve  of  a  common  engine. 


Fig.  2  Ignition  pressure  time  curve. 

In  the  figure,  the  rising  section  can  be  approximately  expressed 
mathematically  as 


Pit) 


When  there  is  only  black  powder  combustion,  the  ignition 
pressure  in  the  combustion  chamber  and  the  production  quantity  of 


the  black  powder  gas  have  the  following  relationship 


_cfm a.  if(o_ 

sc.u  PW+c.l„.  d, 


dP(t) 


After  substituting  in  P(t)  and  — ,  we  obtain 

Cf'  (k)A,  ( _  ZPj\  Z PJ1  \  V,  (_  _6 PJl  jPj_  \ 

m,,~  v  a  ti  )+ctfPi,\  ti  ti  ) 


Because  the  influence  of  environmental  pressure  during 
black  powder  combustion  is  very  small,  the  expression  of  the 
production  rate  is  still  suitable  in  the  combined  combustion 
period. 

(E)  Establishment  of  a  Set  of  Differential  Equations 


By  substituting 

can  obtain 

the  expression  of 

m. 

ig 

into  formula  (2) ,  we 

dP.  C,fk 

\aiFl(k)A, 

(  -2 PJ’ 

3 PJ* ' 

V,  (  6  PJ1 

dt  ~  V , 

*-  V C,{ Pi, 

[  tl  + 

ft-. 

/  c,/P,A  <• 

.  jpj\ + o suP-  C«F w A'P--  + 

+  ti  )+pW‘  —  +  (.c,fky  dt  J 


It  is  now  still  impossible  to  carry  out  numerical  integration 
of  this  differential  equation  because  under  certain  circumstances 
we  only  know  the  function  relation  of  combustion  surface  s  and 
meat  thickness  S'  but  are  unclear  of  the  function  relation  of 
s  to  t.  For  this  reason,  it  is  necessary  to  add  two  sets  of 
equations  to  form  a  set  of  differential  equations  so  as  to  attain 
an  integral  solution.  Because 

db=udt  =  uxPn,dt, 
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therefore 


ds  _  ds  db 
HT~  db  dt 


Integral  initial  values  s  and  (£*  of  the  set  of  differential 
equations  of  the  combined  combustion  period  both  use  the  grain's 
initial  combustion  area  s  and  initial  combustion  meat  thickness 
0=0.  At  present,  the  integrated  pressure  initial  value 
and  the  time  of  ignition  t^  can  only  be  determined  by  an 
approximation  method. 


On  the  basis  of  the  above  given  set  of  differential  equations 
of  the  combined  combustion  period's  pressure  change  pattern,  we 
wrote  a  computer  program  to  predict  the  engine's  initial  pressure 
peak  and  separately  carried  out  prediction  calculations  of  the 
ignition  pressure  peak  of  engines  with  certain  products.  The 
results  and  test  values  are  basically  in  agreement. 

III.  Factors  Influencing  the  Ignition  Pressure  Peak 

It  is  not  difficult  to  see  from  the  above  set  of  differential 
equations  that  in  a  small  solid  rocket  engine  with  very  little 
corrosive  combustion,  the  existence  of  the  production  rate  of 
black  powder  is  the  major  factor  influencing  the  initial  pressure 
peak.  We  call  the  initial  pressure  peak  formed  because  of  this 
the  ignition  pressure  peak  and  call  the  initial  pressure  peak 
formed  from  corrosive  combustion  the  corrosion  pressure  peak. 
Under  normal  conditions,  the  initial  pressure  peak  of  the  engine 
is  composed  of  these  two  parts.  Tests  prove  that  in  small  solid 
rocket  engines  with  rational  designs,  corrosive  combustion  can 
often  be  very  small  and  thus  need  not  be  considered.  However, 
for  a  small  engine  which  uses  black  powder  ignition,  the  ignition 
pressure  peak  that  appears  at  high  temperatures  is  difficult  to 
avoid.  Therefore,  generally  speaking,  the  ignition  pressure  peak 
is  the  main  part  making  up  the  initial  pressure  peak  of  the 
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small  solid  engine  at  high  temperatures. 


This  type  of  pressure  rise  is  usually  expressed  by  peak  value 
ratio  and  at  present  the  high  temperature  peak  value  ratio 

of  most  small  solid  engines  which  use  black  powder  ignition  is  in 
the  1. 4-2.0  range.  This  can  be  even  higher  in  special  situations. 


In  order  to  investigate  the  various  factors  influencing  the 
ignition  pressure  peak  and  find  a  method  for  decreasing  or 
eliminating  the  ignition  pressure  peak,  we  wrote  a  computer 
program  to  carry  out  a  relatively  large  number  of  simulated  cal¬ 
culations  of  the  engine's  ignition  pressure  peak  and  to  study  the 
influence  of  the  engine's  design  parameters  on  the  high  tempera¬ 
ture  ignition  peak.  The  calculation  results  agreed  with  the 
phenomena  observed  in  the  experimental  research  on  the  engines. 


1.  Properties  of  the  Propellant 

When  guaranteeing  conditions  wherein  the  engine  has  reliable 
low  temperature  ignition,  the  larger  combustion  rate  pressure 
index  n  of  the  engine  at  high  temperature,  -the  larger  the  high 
temperature  ignition  pressure  peak.  At  high  and  low  temperatures, 
the  larger  the  difference  of  the  ignition  energy  required  by  the 
engine  on  the  unit  combustion  surface  area,  the  larger  the 
ignition  pressure  peak  of  the  engine  at  high  temperatures. 

However,  because  the  present  propellant  tests  and  employed 
separation  are  more  serious,  this  problem  has  yet  to  receive 
enough  attention  and  it  is  still  not  regarded  as  an  important 
performance  indication  of  the  engine.  This  is  worthy  of  attention. 


2.  The  Loading  Density  of  the  Engine 

When  guaranteeing  conditions  of  reliable  low  temperature 
ignition  in  engines  which  use  the  same  solid  propellant,  the 
larger  the  loading  density  of  the  propellant  in  the  combustion 
chamber  (the  smaller  the  initial  free  volume) ,  the  higher  the 
ignition  pressure  peak  of  the  engine  at  high  temperatures. 
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Therefore,  as  regards  a  certain  type  of  propellant,  in  order  to 
prevent  the  formation  of  an  excessively  large  ignition  pressure 
peak  at  high  temperatures,  we  should  limit  the  loading  density 
of  the  propellant  in  the  combustion  chamber. 

In  order  to  resolve  the  problem  of  a  certain  engine  having 
a  high  initial  pressure  peak  at  high  temperatures  and  sunken 
initial  section  pressure  at  low  temperatures,  we  carried  out  a 
large  number  of  ignition  tests.  The  results  proved  that  at  high 
temperatures  the  engine's  initial  pressure  peak  noticeably 
changed  with  the  loading  density  of  the  combustion  chamber.  At 
the  same  time,  the  test  phenomena  also  showed  that  after  the 
initial  free  volume  of  the  engine  was  smaller  than  a  certain 
determined  value,  at  high  temperatures  the  engine's  initial 
pressure  peak  drastically  changed  with  the  reduction  of  the 
initial  free  volume.  When  the  initial  free  volume  changed  in 
another  relatively  large  range,  at  high  temperatures  the  engine's 
initial  pressure  peak  changes  were  relatively  mild.  Therefore,  in 
engine's  with  relatively  large  loading  densities,  we  should 
control  the  tolerance  of  the  combustion  chamber's  initial  free 
volume  and  then  we  will  be  able  to  ensure  that  the  engine's 
initial  pressure  peak  is  stable  at  high  temperatures. 

Although  the  loading  densities  of  propellants  in  greatly 
reduced  engine  combustion  chambers  can  eliminate  the  ignition 
pressure  peak  at  high  temperatures,  yet  this  will  cause  the 
loading  quantity  of  the  propellant  to  greatly  decrease.  Therefore, 
when  using  this  method  to  decrease  the  engine's  ignition  pressure 
peak,  we  must  consider  it  comprehensively. 

3.  Design  Parameter  of  the  Engine  -  Surface-to  Throat  Ratio 
Rot 

It  is  common  that  when  the  other  parameters  of  the  engine  are 

constant,  the  larger  the  surf ace-to- throat  ratio  R  ^ ,  the  smaller 

ot 

the  peak  value  ratio  at  high  temperatures.  However,  an  excessively 
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large  R  causes  the  maximum  pressure  to  increase  and  thus 
increases  the  negative  weight;  on  the  contrary,  if  the  R 
excessively  small,  this  can  only  cause  V*  to  markedly  increase 
but  can  also  influence  the  release  of  propellant  energy  because 
the  equilibrium  pressure  is  excessively  low.  Therefore,  as 
regards  each  type  of  solid  propellant,  when  designing  engines, 
it  is  necessary  to  select  a  moderate  R  value  under  totally 
allowable  conditions  and  this  will  then  be  able  to  benefit  the 
raising  of  the  overall  performance  of  the  engine. 

4.  Proper  Ignition  Design  Has  Certain  Effects  on  Lowering 
the  Engine's  High  Temperature  Ignition  Pressure  Peak 

(A)  Use  of  High  Energy  Ignition  Powder 

High  energy  ignition  powder  generally  has  a  low  combustion 
rate  pressure  index,  small  gas  content  and  high  heat  value,  it 
has  very  strong  energy  output  in  an  infrared  area,  it  can  have 
various  production  prescriptions  and  grain  compositions,  it  is 
an  effective  igniter  and  there  is  a  tendency  under  many  conditions 
to  replace  the  black  powder  with  a  high  energy  igniter,  its  use 
can  to  a  certain  extent  lower  the  engine's  ignition  pressure 
peak  and  improve  the  low  temperature  ignition  performance. 

(B)  Revising  the  Structural  Form  of  the  Ignition  Design 

The  ignition  energy  required  on  the  ignition  unit  surface 

area  for  different  types  of  propellants  under  the  same  initial 
temperature  are  different.  At  the  same  time,  during  the  ignition 
process,  the  heat  response  time  required  to  establish  a  normal 
heating  layer  on  the  propellant's  surface  is  also  different  (for 
the  majority  of  propellants,  the  suitable  response  time  is  usually 
in  the  5-20ms  range) .  When  the  ignition  powder  is  a  certain 
quantity  and  the  heating  time  of  the  black  powder  gas  is  either 
too  short  or  too  long  on  the  propellant's  surface,  these  are 
disadvantageous  to  the  normal  ignition  on  the  propellant's 
surface.  When  carrying  out  ignition  design  of  engines  with 
different  types  of  grains,  we  should  distinguish  the  different 


conditions  and  suitably  select  the  black  powder's  quantity,  size, 
type  of  igniter  and  location  of  the  igniter  etc.  so  as  to  obtain 
optimal  ignition  effects. 

If  static  tests  of  a  certain  engine  show  that,  given  the 
other  conditions  are  the  same,  when  using  a  celluloid  ignition 
box,  the  bursting  pressure  of  the  box  is  low  and  the  initial 
section  of  the  engine's  P-t  curve  has  a  relatively  smooth  rise. 
When  using  an  aluminum  ignition  box,  the  bursting  pressure  of 
the  box  is  high,  the  initial  section  of  the  engine's  P-t  curve 
assumes  a  saddle  shape  and  the  maximum  pressure  at  high  tempera¬ 
tures  is  about  10-15%  higher  than  when  using  a  celluloid  box. 

We  can  see  that  good  ignition  design  has  certain  effects  on 
lowering  the  ignition  pressure  peak  of  engines  at  high  tempera¬ 
tures. 

IV.  Experimental  Research  on  Eliminating  the  Ignition  Pressure 
Peak 

On  the  basis  of  specialized  experimental  research,  we  propose 
a  new  ignition  design  method.  In  small  solid  engines  with  not 
very  serious  corrosive  combustion,  use  of  this  method  can 
effectively  eliminate  the  initial  pressure  peak  of  the  engine  at 
high  temperatures  under  conditions  wherein  the  engine  grain  is 
not  reduced  and  there  is  no  increase  of  the  negative  weight  and 
thus  improve  the  overall  performance  of  the  engine. 

The  procedure  of  the  new  design  method  is  as  follows: 

1.  In  developing  the  grain  engine,  we  begin  from  the  require¬ 
ments  of  low  temperature  reliable  ignition  and  not  having  a 
relatively  large  ignition  pressure  peak  at  high  temperatures,  and 
carry  out  predictions  and’  calculations  of  the  new  ignition 
material  and  structural  design. 

2.  When  using  false  grain  in  the  engine,  we  measure  the 
actual  ignition  pressure's  time  curve  in  the  engine  to  examine 
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whether  or  not  it  has  reached  the  predicted  design  requirements 
of  the  original  igniter.  It  is  necessary  to  make  revisions  on 
the  original  ignition  design  until  the  ignition  design  reaches 
design  requirements. 

3.  On  the  basis  of  the  two  operations  mentioned  above,  we 
carry  out  ignition  tests  of  the  real  engine  and  measure  the 
P-t  curve  of  the  engine. 

This  new  ignition  design  method  has  already  been  applied  in 
certain  engines.  Figures  3  and  4  separately  show  the  P-t  curves 
of  engines  at  high  and  low  temperatures  when  using  the  new 
ignition  design  and  when  using  the  common  ignition  design.  It  is 
not  difficult  to  see  from  the  curves  that  in  using  two  different 
types  of  ignition  designs  in  the  same  engine,  when  we  are  able  to 
guarantee  low  temperature  reliabble  ignition  and  stable  combus¬ 
tion,  there  is  a  marked  difference  in  the  pressure  curves  of  the 
engine  at  high  temperatures.  When  using  the  new  ignition  design 
method,  the  peak  value  of  the  engine's  pressure  curve  at  high 

temperatures  is  relatively  small,  Y"  =1 . 08 (approximately) ,  while 

P 

when  using  the  common  ignition  design,  the  peak  value  of  the 
pressure  curve  at  high  temperatures  is  more  than  1.56  higher 

than  T  . 
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Fig.  3  Pressure  time  curves  of  a  certain  engine  when  using  the 
new  ignition  design. 
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Fig.  4  Pressure  time  curves  of  a  certain  engine  when  using  the 
common  ignition  design. 

Comparisons  of  ignition  tests  on  real  engines  show  that  the 
new  ignition  design  method  has  marked  effects  on  raising  the 
overall  performance  of  engines,  it  has  a  simple  structure, 
reliable  operations  and  use  is  convenient  etc.  It  can  be  used 
on  small  solid  engines  and  it  is  also  easy  to  realize  on  rocket 
weapons  with  finalized  designs. 
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